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»"~~ | Abstract. A general dynamically consistent 2D flat distribution function is built to model the kinematics of neigh- 

[j bouring stars. Application leads to the measurement of a short galactic scale length R p between 1.7 and 2.9 kpc and 

a locally decreasing rotation curve. This is in agreement with recent determinations based on kinematics and counts 

of distant stars. These results rule out the classical assumption that 2R p = R a or that a z (R)/an(R) is constant when 

the vertical scale height h z (R) is assumed to be constant. We explain why the measured squared axis ratio of the 

J> ■ velocity dispersions cr^/cr^ of disc stars is less than 1/2. This ratio has been claimed to be important evidence for the 

CO | non-axisymmetry of the galactic disc. We show that this is not the case and that it may be simply explained with 
a realistic axisymmetric disc model if the circular velocity is locally declining or if there is a mismatch between the 

^j \ photometric and kinematic scale lengths. 
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^ ' 1. Introduction 



Dynamical self-consistent constraints deduced from the Boltzmann equation are necessary to analyse and interpret 
kinematic data in the Galaxy. It has been shown that the standard epicycle theory approach, a first order theory to 
estimate galactic characteristics, is subject to systematic errors (see Kuijken & Tremaine 1991, who put forward higher 
order developments). In some recent work, distribution function solutions of the Boltzmann equation (proposed by 
Shu 1969) and/or third moments in Jeans equations (Cuddeford & Amendt 1991) have been used to obtain a more 
consistent description of our Galaxy (Cuddeford & Binney 1994, Kuijken & Tremaine 1991, Evans & Collett 1993, Fux 
& Martinet 1994). With different approaches, they attempt to determine and to make use of exact constraints between 
global parameters describing our Galaxy, namely the slope of the circular velocity curve, the density and kinematic 
scale lengths etc... 

In this paper, we give a short review of observational constraints and the different values taken by these authors. 
Secondly we build a new model of Shu-type distribution functions where scale lengths and the shape of the velocity 
curve are free parameters. Thirdly we compare and adjust this model to velocity distributions of stars in the solar 
neighbourhood taken from the Gliese catalogue and other catalogues. Implicit constraints are related to the asymmetric 
drift relation and the Lindblad equation that gives the ratio of radial and tangential velocity dispersions. We conclude 
that extended Shu distribution functions recover most of the fundamental kinematic properties of the galactic disc. We 
show that it is possible to build a consistent model explaining local stellar kinematics and to determine local structural 
galactic parameters. 
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2. Observational constraints 

We use local kinematic data from Gliese & Jahreiss (1991) and other catalogues to constrain galactic structure. Stars 
near the sun with large relative velocity provide information on non-local galactic structure like the kinematic and 
density gradients or the slope of the rotation curve. The local star velocity distribution is a section of phase space 
and this observed distribution depends on non-local quantities like the potential, or global density distribution of 
stars (other examples are given by Dehnen and Binney 1995, who describe ways to map the Galaxy's gravitational 
potential and distribution of matter). We show that the analysis of local kinematics is able to recover the potential 
gradient and stellar scale lengths. In previous analyses, assumptions have been either a flat velocity curve (Evans & 
Collett 1993) or a relation between density and kinematic gradients (Cuddeford & Binney 1994). These assumptions 
allow these authors to obtain tractable distribution functions. However they found that these distribution functions 
are not compatible with existing observational evidence. We discuss observational evidence on galactic structure, and 
emphasize the existing systematic differences between authors. 

2.1. The circular velocity curve 

Determinations of the galactic circular velocity curve are based on radial velocity observations of objects belonging 
to the disk population. Due to their small velocity dispersion, they follow the galactic rotation curve closely. However 
exact determination of the velocity curve depends critically on the adopted galactic radius Rq and circular velocity <d 
at the solar position. Assuming Ro = 7.9 kpc and 9o = 185 km s _1 , Rohlfs and Krcitschmann (1988) found that the 
velocity curve "is slowly declining from 9 = 200 km s^ 1 at R = 4.5 kpc to 9 = 173 km s" 1 at R = 11 kpc" (however 
taking 9o = 220km s _1 , they obtained a flat rotation curve). This result is also corroborated by Fich et al. (1989), 
Pont ct al. (1994) and by Dambis et al. (1995). 

Since most recent determinations favor values around Rq = 7.9 kpc (Reid 1993) or Ro — 8.1 kpc (Pont et al. 1994) 
and 9o = 185 km s , much smaller than the recommended values of Rq = 8.5 kpc and 9o = 220 km s _1 by the 
IAU general assembly in 1985 (Kerr & Lynden-Bell 1986), we consider it very likely that the circular velocity curve 
is locally declining. We estimate a ~ —0.1 to — 0.3 for v c (R) = R a at the solar galactic radius of Ro = 7.5 kpc and 
9 = 185 km s" 1 . 

Considering the Lindblad equation ai/a^ ~ 2 (1 H — dinR ) ' we concm de that there is no a priori conflict between 
the slope of the rotation curve and the observed ratio of velocity dispersion o~\I<j\ smaller than ^ (Kerr & Lynden-Bell 
1986). 

2.2. The disc density scale length 

Most disc density scale length determinations range between 2.5 and 5 kpc (see Kent et al. 1991 and Robin ct al. 1992a), 
but part of these measurements are model dependent and require assumptions. Assuming uji(R) = exp(— R/2h), Lewis 
& Freeman (1989) gave h = 4.4 kpc from velocity dispersions of distant K giants, van der Kruit (1986) obtained a 
larger value of 5.5 kpc from the Pioneer 10 background experiment. But his determination gives access only to the 
radial to vertical scale ratio, and he assumed a h z — 325 pc vertical scale height for the old disc. Adopting the recent 
determination of 250 pc for h z (Kuijken & Gilmore 1989, Haywood et al. 1996), the radial scale length deduced from 
Pioneer 10 data would be closer to 4.2 kpc. Recent direct determinations range between 2.5 kpc (star counts: Robin 
et al. 1992ab, Ojha et al. 1996, COBE map: Durand ct al. 1996) and about 3.5 kpc (Kent et al. 1991). We note that 
the kinematic determination of the radial scale length obtained by Fux & Martinet (1994) is h — 2.5 kpc, assuming 
a radially constant scale height h z , and h = 3.1 kpc with a positive local h z gradient of 30pc/kpc. Their analysis is 
based on the asymmetric drift equation and on moments of the Boltzmann equation. 

2.3. The disc kinematic scale length 

We know only few direct determinations for the kinematic scale length of stellar discs in the Milky Way. Neese & Yoss 
(1988) measured the gradient of radial velocity dispersion from 364 stars mainly towards the galactic anticentre and 
found da R /dR = -3.8 ± 0.6km s" 1 kpc" 1 or d\na 2 R /dR = -(SJkpc)" 1 . Lewis & Freeman (1989) obtained 4.4kpc 
for the o\ scale length (or 8.8 kpc for o-r). These results are based on 600 distant giants towards the galactic centre 
and anticentre. As remarked by Evans & Collett (1993), the Lewis & Freeman (1989) fit is probably not good at the 
solar position and a different value for the velocity dispersion at the solar galactic radius based on K giants (Delhaye 
1965) should be used. With this value Evans & Collett (1993) tried a non-exponential function to model the radial 
dependence of velocity dispersions and obtained a better fit. Using this function, we obtain for the local kinematic 
scale length dh\a\/dR = (— 4.8kpc) _1 = — 0.21kpc _1 . This illustrates the precision of the determination. From 
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star counts and proper motions Ojha ct al. (1996) measured velocity dispersion and gradients and obtained similar 
kinematic scale lengths. 

It is usually accepted that galactic discs have constant scale heights and exponentially decreasing velocity dis- 
persions. Measured velocity dispersions in galactic discs led Bottema (1993) to conclude that the velocity dispersion 
decreases exponentially with radius like a\ ~ e~ R / Rp if the vertical to radial velocity dispersion ratio is constant 
with radius. This last current assumption (or the plane-parallel assumption) must be ruled out if, for example, the 
observations in the Milky Way show that R p ^ R a i . 

For a given disc galaxy, there is not a strong chromatic dependence of the scale length (Elmegrccn & Elmegrecn 
1984, Giovanardi & Hunt 1988). The bulk of 86 spiral galaxies (Fig. 9 by dc Jong, 1996) shows R b /Rk ~ 1-2. This is 
a strong indication that the stellar populations are distributed in discs with similar scale lengths. For our model fitting 
in the next section, we will assume no chromatic dependence for our Galaxy. We will assume the same for kinematic 
scale lengths though there are no existing observations that confirm this. 

2.4- Disc axisymmetry 

Evidence from a variety of sources shows the existence of a bar in the central few kpes of the galactic center, with 
rough agreement on the orientation of the bar. However it is much more difficult to find evidence from star counts 
in the galactic plane for ellipticity of the stellar disc at the solar galactic radius, particularly if we are on one of the 
symmetry axes. 

Signatures from kinematics give more reliable constraints (Kuijken & Tremaine 1994): 1) Oort's constants C and 
K and the LSR radial velocity show no evidence for local non-axisymmetry, 2) while vertex deviation for low velocity 
dispersion and young stars may show such evidence this last effect could be caused by very local non-stationary effects 
related to the presence of spiral arms, 3) for stars in the solar neighbourhood of high velocity dispersion, the velocity 
ellipsoid points towards the galactic centre (I — 5.5° ± 4.2), with no indication of local non-axisymmetry, 4) the axis 
ratio <J V /<T U is observed to be about ~ 0.5 and can be a signature of ellipticity but as we show in detail in this paper, 
it may also be explained if the rotation curve is locally decreasing or if the kinematic and density scale lengths are 
very different. 

Stronger arguments in favour of the effect of a bar at R = Rq come from the Blitz & Spergel (1991) analysis of 
HI data from which they deduced that the LSR has a radial motion of 14 km s _1 . However the effect of the warp 
beyond R — 12kpc (Burton & te Lintel Hekkert 1986) may modify their analysis and the LSR radial motion is not 
confirmed by other observational measurements: v r (LSR) = — 1 ± 9 km s _1 (Table 3 in Kuijken & Tremaine 1994). 
As a conservative hypothesis we will still assume that the galactic disc is axisymmetric and we will show that local 
kinematics may be explained in the frame of this hypothesis. 

2.5. The vertical structure 

A plane-parallel potential is developed in the next section. Such potentials do not permit the consideration of the 
correlations between vertical and radial motions that must exist. From the probable mass galactic distribution, these 
correlations have certainly a better representation with a spherical potential than with a plane-parallel one (this is not 
true for the density distribution). Such a spherical potential will be briefly developed in the next section and general 
solutions fitting observational constraints will be given in Section 4. Intermediate galactic potentials between plane 
parallel and spherical do not have 3 known integrals of motion (with the exception of Stackel potentials) and cannot be 
used to find simple analytic distribution functions. For this reason we restrict the present analysis to Stackel potentials 
to include the effect of the vertical structure. In these cases, the distribution functions are deduced following Statler 
(1989) and are extensions of the Shu and Schwarzschild distribution functions. 

Vertical structure introduces a bias on the Lindblad and asymmetric drift equations which is estimated by Cud- 
deford & Binney (1994). The complete analysis done by Fux and Martinet (1994) allows them to measure locally the 
third order derivative term of the potential related to the flatness of the potential. This derivative is related to the 
radial gradient of the vertical scale height h z . 

3. Galactic disc distribution function 

Shu (1969) described a phenomenological distribution function (DF) that has the great advantage of relative simplicity 
and that has been used frequently. This distribution function is a two dimensional representation of a flat and rapidly 
rotating stellar disc in an axisymmetric potential (rotational symmetry). The DF is written explicitly in terms of 
two integrals of motion (energy and angular momentum) and is a steady-state solution of the Boltzmann equation 
if the potential is time- independent. For small velocity dispersions, Shu's DF closely approximates the Schwarzschild 
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velocity ellipsoid DF. Two functions are free: the input surface density distribution S(-R) and the input radial velocity 
dispersion a(R) (Shu 1969, Tamisier 1991). The main advantage of the Shu DF is that the associated exact surface 
density distribution S(i?) and radial velocity dispersion o~r(R) (DF moments of and 2 nd order) will in general remain 
close to the input functions S(i?) and a(R), even with quite large velocity dispersions. We define the Shu DF as 

,, p r , 2Cl(R c ) S (-R c ) r (E — E arc (R c )) 
nE ' L ' )= k(R c )2^{R c ) CM ^ (R c ) ] lf ^ >0 (1) 

f(E,L z ) = if L z <0, (2) 

where R c = R C (L Z ) is the radius of the circular orbit with angular momentum L Zl (R C (L Z ) is bijective for potentials 
considered in the next sections), E C i rc (R c ) is the energy of a circular orbiting star at radius R c with angular momentum 
L z . il, is the angular velocity and k is the epicyclic frequency 2f2[l + |d InQ/d In R] 1 / 2 . The exponential part of the 
DF (Eq.l) represents a ring of stars with angular momentum L z rotating near the radius R c and with radial velocity 
dispersion a. The other terms are normalizations which permit a distribution of elementary rings with different L z 
(and R c ) and with a total mean density close to S(i?). Equation 2 is given to ensure that the total angular momentum 
of the disc is non-zero. In this model the disc has a maximum rotation compared to any other distribution function 
with the same resulting density and no retrograde orbits are allowed. Here these restrictions have no consequence since 
we analyse the local kinematics where only one retrograde orbit is found among one thousand stars. 
Since we want to build a model whose resulting density is S(i?) ~ So cxp(— R/R p ), we put 

Z(R C {L Z )) = S exp[-Rc(L x )/R p ] (3) 

and in order to have a nearly exponential o~r{R) distribution, we put 

a(R c {L z )) = a cxp[-R c (L z )/R a }. (4) 

Thus the DF depends only on the integrals of motion and on the following free quantities: the velocity curve v c (R) 
and the two constant scale lengths R p and R a . 

3.1. Flat rotation curve model 

We describe a model with a flat rotation curve v c (R) = 1 and a potential given by &(R) = ln(i?). Then we have 
Rc(L z ) — L z /v c (R c ) = L z . Pertinent parameters are the amplitude (To and the scale length ratio R a /R p . We define 
models with quasi-exponential density and dispersion distributions: 

Z(R C {L Z )) = S exp[- L z /R p ] (5) 

and 

a(R c (L z )) = o-o exp[-L z /R a }. (6) 

For different values of cro, Figs. la-2a show that the density £(i?) and the radial velocity dispersion cfr{R) remain 
close to an exponential (with scale lengths R p and R a respectively) for models with the smallest velocity dispersions. For 
the largest dispersion (<7o = -64), the effective density scale length is larger than the R p parameter. The mean streaming 
velocity v<p (Figs. lb-2b) of the modelled stellar populations departs from the circular velocity curve (v c (R) = 1). This 
is due to asymmetric drift which is directly related to the amplitude of the radial velocity dispersion. For the models 
with large cro, the ratio o^Jgr (Figs. lc-2c) departs from the limiting value \[2j2 for zero velocity dispersion. 

These results can be directly compared to that of Evans & Collett (1993) and Cuddeford & Binncy (1992): the 
velocity dispersion ratio can be much smaller than -\/2/2, as observed in the solar neighbourhood (Table 8 of Kerr & 
Lyndcn-Bell 1986) in at least three situations: Ro is smaller than R p , R a is much larger than R p or the rotation curve 
is decreasing with radius. 

For small velocity dispersions, the ratio o~(j,/o~r depends only on the shape of the velocity curve. This is not true 
for large dispersions, where dependence on the scale lengths is significant. 

Shu-type models exist with the Lewis & Freeman (1989) assumption that 2R p = R a (Fig. 1), (see for example 
Binncy 1987, Kuijken & Tremaine 1991 and Cuddeford & Binney 1992). Such models imply that R p > 5kpc (i.e. 
Ro/Rp < 1.6) to explain measured ratios o-^/ctr smaller than \[2j2. It may be also accounted for if R„j ' R p >> 2 (Fig. 
2)- 
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Fig. 1. (a) The surface density E(-R) and the radial velocity dispersion <jr(R) generated in a logarithmic potential by the 
distribution of Eqs. 1, 7 and 8 with R a = 2R P and four different cro =(.08, .16, .32, .64). (b) The mean streaming velocity, (c) The 
velocity dispersion ratio. 
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Fig. 2. Same as Figure 1, but with R a — 6R P . (a) Surface density and radial velocity dispersion, (b) The mean streaming 
velocity, (c) The velocity dispersion ratio. 



3.2. Power law rotation curve models 

In what follows we use a more general model with potentials yielding a power law circular velocity curve: v c (R) = R a . 

1 

The radius of a circular orbit with momentum L z is given by R C (L Z ) = L z /v c = L z +1 and the circular velocity is 
v c {L z ) = L z +1 . This gives for S and a: 



E(L Z ) = S exp[-L° +1 /-Rp] 



a(L z ) = a exp[-L° +1 /i? CT ] 



(7) 
(8) 



Here again the resulting density and dispersion will be in general nearly exponential with scale lengths R p and R a , 
but models loose this convenient behaviour when the velocity dispersions are too high. 

Other Shu models could be built using different closed forms or using some spline functions to fit data obtained 
along the galactic radius. Such phcnomcnological models are certainly less elegant than for example Evans & Collctt 
(1993) exact models that give DFs for a few defined density or dispersion distributions. However Shu models arc more 
flexible, cover a larger range of model types and they have a sufficient number of free parameters for basic comparison 
to observations. Thus they allow a more extensive comparison to available data. 
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3.3. 3D models 

3D models arc built to include the vertical structure of the potential. We detail two cases for which explicit DFs are 
easy to obtain. In the case of a plane parallel potential, the DF may be multiplied by the trivial term: 

j2^^T7) exf ^ ' (9) 

For a spherical potential, a 3 dimensional DF may be deduced from the DF given by Equ. (1) and from the equation 
lib given by Statler (1989) for Stackel potentials; in the case of a fiat rotation curve, the vertical term is: 



(2tt)i/2S;(L 2 



-exp 



-+(: 



(10) 



This shows that the correlation term is the most different for the highest velocity dispersion populations. In the next 
section we consider plane parallel and spherical models with power law rotation curves. Other more general models 
using DFs in Stackel potentials (Statler 1989) have also been used to examine solutions in intermediate potentials 
between plane-parallel and spherical ones. In the present context DFs in Stackel potentials introduce a new parameter: 
the focus of the ellipsoidal system of coordinates (Statler, 1989). For all these models we have assumed that the (radial) 
vertical and (radial) horizontal dispersion kinematic scale lengths are equal R„ % = R a . The results of these intermediate 
models are presented in Section 4.3 and can be compared to the plane-paralel model solutions. 

4. Models versus Observations 

4-1. Observational data 

In this section we compare the models presented in the previous section to the local stellar kinematics as given by two 
catalogues. The first one is the recently published Reid et al. (1995) catalogue that improves considerably the data for 
K and M stars found in the CNS3 catalogue of nearby stars (Gliese & Jahreiss 1991). They include new radial velocity 
measurements and reject stars with incorrect spectral types or giants previously misclassified as dwarfs. 

We build a second catalogue by extending as much as possible the CNS3 catalogue using the Simbad database at 
the CDS in Strasbourg and by completing missing data, colours, radial velocities and some trigonometric parallaxes 
when new observations are available (mainly from the Hipparcos Input Catalogue, Turon et al. 1992) but also from 
a bibliography of radial velocities (Barbier-Brossat et al. 1994) and the General Catalogue of Trigonometric Stellar 
Parallaxes (Van Altena et al. 1991). We remove misclassified stars identified by Reid and stars with M v > 8. We retain 
only O to G type stars to get a catalogue complementary to Reid's catalogue. 

We extract from these catalogues distance-limited samples in order to obtain homogeneous and kincmatically 
unbiased 3D velocity distributions. We use distance criteria given by Reid et al. (1995) for their catalogue, and apply 
a similar process to define the completeness limits for the second catalogue (26 pc for M v < 3, 23 pc for 3 < M v < 5 
and 21 pc for 5 < M v < 8). 

Finally we consider multiple stars as one kinematic object since the kinematics of the members are correlated. We 
reject one star with a velocity modulus larger than 200 km s _1 that is certainly a halo star unmodelled by our DF. 

Stars or data not in the CNS3 and Reid catalogues are given in Tables 1 and 2. 

4-2. Model adjustment 

Distribution functions are modeled by Eqs. 1, 2, 7 and 8. Data-model comparison is done by adjustment of model 
parameters using the maximum likelihood method in the (vr, Vcf,) space for plane-parallel models (and vr, v^, v z 
for 3D models). All stars are assumed to be at solar galactic radius R = Rq. The Uq, vq components of the solar 
velocity relative to the LSR are two local parameters. Other parameters S , oo are dependent upon the local number 
of stars and the local radial velocity dispersions. Non-local parameters arc the kinematic and the density scale lengths 
of stellar populations and the slope of the circular velocity curve. These quantities are obtained as a mean over a few 
kiloparsecs around the solar galactic radius corresponding to the extent of the radial excursions of stars that are in the 
solar neighbourhood, and that collectively carry some information from these more distant regions. For 3D models, 
the coordinate of the focus of the ellipsoidal coordinate system used to define the Stackel potential is a supplementary 
adjustable parameter. 

The radial to tangential velocity dispersion ratio of a low dispersion population constrains mainly the slope of the 
circular velocity curve. For kincmatically hot populations, this ratio depends also on R p and R a - The mean velocity 
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Table 1. Stars with new measurements labelled with *. Other data are taken from the CNS3. New V ra d are extracted from 
Barbier-Brossat et al. (1994), trigonometric parallaxes from the General Catalogue of Trigonometric Stellar Parallaxes (Van 
Altena et al. 1991) and B-V from the Hipparcos Input Catalog (Turon et al. 1992). 
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Table 2. Supplementary stars. Data obtained or deduced from the Hipparcos Input Catalog (Turon et al. 1992). 
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Utj, is determined by asymmetric drift. Comparing stellar groups with small and large velocity dispersions gives a more 
accurate measurement of asymmetric drift. For this reason we divide the catalogue in samples (identified by a number 
k) according to the absolute magnitudes. The sample including the smallest absolute magnitude corresponds to the 
brightest stars or the bluest since all these stars are on the main sequence. This sample includes young stars and 
has their kinematic characteristic: low velocity dispersion. In contrast, samples with fainter stars have redder main 
sequence stars and are a mixture which includes older stars with higher velocity dispersion. We determine the sample 
division by minimizing the errors on the R p and R a estimates. 

Maximum likelihood is used to estimate the best model parameters. The likelihood function L is the product of 
likelihoods L&, where Lk is the likelihood of sample number k. 



Lk = II fk( v R,i, v 4>,i)/^k 
i=l,N k 

with 

Sfe = / / fkdv R dv<f, 

where N^ is the number of stars in sample k. 



(11) 



(12) 
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Fig. 3. Radial and tangential velocity histograms for three samples of stars sorted according to their absolute magnitude M v . 
The maximum likelihood models are plotted (continuous lines) as well as individual DFs (dashed lines). 

Each sample is modeled by using two elementary distribution functions, fk = fk,i + fk.2 (continuous lines in Fig. 
3) corresponding to the sum of two DFs (dashed lines in Fig. 3) given by Eqs. (1, 2, 7 and 8). 

Jf.,3. Results 

Solutions do not depend significantly on the catalogues, the splitting method or the shape of the potentials considered 
(spherical, plane-parallel or intermediate) . Best fits with the smallest errors arc obtained using both catalogues simul- 
taneously. Figure 3 shows the observed velocity histograms and the corresponding model curves for the best estimates 
of parameters with a plane-parallel potential and i?o = 8.5 kpc and Vo = 220 km s _1 . This model is in good agreement 
with data (the reduced x 2 is approximately one for each histogram). 

The results obtained with each catalogue are given in Table 3 in case of plane-parallel models and in case of interme- 
diate 3D models. R p , e ff and R a ,eff scale lengths are the effective scale lengths of the computed density or dispersion 
that, in the case of high velocity dispersions, differ from the R p or R a parameters used in Eqs. 3 and 4. Formal errors 
are small but we find a strong correlation of scale length with solar velocity vq which limits the accuracy of density 
scale length estimates and could bias our conclusions. 

4-4 ■ Sample splitting 

We divided the stellar catalogue into one, two or three samples with various sizes and we made systematic explorations. 
We found subdivisions that minimize errors of non-local parameters. Within the error bars, other subdivisions gave 
solutions in agreement with the best fit. The best fit is obtained with two small samples with a low velocity dispersion 
and a large one with a high velocity dispersion. This corresponds certainly to the best way to constrain the asymmetric 
drift. 



4-5. R p and the solar velocity Vq 

The maximum likelihood solution gives a small value of R p , e ff — 1.7 kpc. Plotting the asymmetric drift relation v'q = 
Vq + c < II 2 > with the CNS3 catalogue yields a result different from Delhaye (1965) who obtained vq = 12 km s _1 . 
It appears that the mean velocity v of the few tens of bluest CNS3 main sequence stars is small and favours a small 
Vq value between 3 and 7km s . This explains the small R p found which is strongly correlated with Vq. Most recent 
vq determinations are small, Gomez and Mennessier (1977) found 6 km s -1 , Mayor (1974) obtained 6.3kms~ 1 and 
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Table 3. Model parameters: solar velocities, density and kinematic scale lengths, slope of the galactic velocity curve. Galactic 
structural parameters are obtained for various splittings of the data in absolute magnitude and assuming Ro — 8.5 kpc and 
V = 220 km s _1 . 

Catalogue number splitting mq Vq Rp,eff R<?,eff a 

of stars km s _1 km s _1 kpc kpc 
Plane-parallel models 

Reid et al. 477 I 8.4 ± 1.5 4.3 ± 4.4 L9 36 -.21 ± .07 

Blue nearby stars 504 2(228/276) 8.7 ±1.1 2.5 ± 2.1 1.8 24-30 -.23 ±.06 

Reid and Blue stars 981 3(228/256/497) 8.6 ± 0.9 2.2 ± 1.7 1.7 ±.3 > 30 -.22 ±.05 

Reid and Blue stars with 981 3(228/256/497) 8.6 ± 0.9 2.6 ± 1.6 2.0 ± .3 > 30 -.18 ±.04 
Ro = 7.9 kpc, Vb = 185 km s" 1 



Reid and Blue stars : 


: Vq 
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8.7 ±.9 
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2.9 ±.3 


18-30 


-.20 ±.05 


3D models 


Reid and Blue stars 
Reid and Blue stars : 


: Vq 


= 8 


981 
981 


3 (228/256/497) 
3 (228/256/497) 


8.5 ±1.2 

8.6 ±1.2 


2.5 ±2.2 
8 


1.7 ±.3 

2. 8 ±.4 


> 30 
20-30 


-.22 ±.06 
-.20 ±.06 



Delhaye (1982) found 8.5km s _1 . Some other references may be found in Kuijken and Tremaine (1991). Oblak (1983) 
built new samples carefully determining stellar ages in order to analyse the asymmetric drift. He found a small solar 
velocity v & — 5.0 ± 0.7 km s _1 . 

To summarize, if we admit 2 and 8 km s _1 as the lower and upper limits on the w© solar velocity, we obtain limits 
on the effective density scale lengths R P , e ff from 1.7kpc to 2.9 kpc. 

The scale length we have determined is proportional to the assumed solar galactic radius Ro', with our best fit it 
gives R p = 1.7 * (Ro/8.5) kpc. 

5. Discussion and conclusion 

Dynamical analysis of the stellar kinematics in our Galaxy is frequently based on the comparison between observed 
and modelled moments of the velocity and not directly on the measured distribution function. A practical reason for 
that is the apparent lack of simple dynamical models predicting realistic distribution functions. Moments contain most 
of the dynamical information. However their measurements based on observational data may be strongly biased by a 
few stars (high velocity members in a binary or halo stars when analyzing the disc). For example Cuddcford & Binncy 

(1992) show how these biases can be reduced. Evans & Collctt (1993) obtain new realistic DFs by generalising the 
Rybicki disc model, but in the case of exponential disk models with R p = R a i , they only obtain moments. 

Here we have shown that the Shu (1969) models that are exact solutions of the collisionless Boltzmann equation, 
can be used successfully to analyse galactic dynamics. We build a nearly exponential stellar disc (much more accurately 
exponential than galactic discs really are), which includes the basic necessary parameters, the density and kinematic 
scale lengths, the slope of the rotation curve and which removes all restraining assumptions done in previous work. 
The comparison to the velocity distribution from an unbiased kinematic sample of nearby stars gives an accurate 
determination of the galactic structure: assuming Ro = 8.5 kpc and Vq — 220 km s _1 , we conclude that the stellar 
density scale length of the Milky Way ranges between 1.7 and 2.9 kpc. We obtain a large kinematic scale length; 
this is in poor agreement with the galactic R a determination obtained by Lewis and Freeman (1989) or by Bottema 

(1993) in external similar spiral galaxies. The large kinematic scale length that we find may result from the large 
proportion of young stars in our sample with a kinematic behaviour resembling gas that has a null kinematic gradient. 
It may result from the inadequacy of the model to also measure this parameter, e.g. if the accuracy of velocities in the 
Gliese catalogue is not sufficient. This could be checked swiftly with Hipparcos data. We find that the rotation curve 
v c (R) — R a is decreasing at the solar galactic radius with a — — 0.22± .05. All these quantities are obtained from local 
stars that have non-negligible radial motions. It means that these structural parameters are mean values measured 
over an extent of a few kiloparsecs around the Sun. Mayor & Oblak (1985) and Oblak & Mayor (1987) proceed partly 
with a similar approach (maximum likelihood and Monte-Carlo) but based on numerical integration of stellar orbits 
in different galactic potentials. They find R a i = 5 kpc or R a = 10 kpc. 

Few other previous kinematic determinations of the Galaxy scale lengths have been obtained assuming R p = R a 2 . 
Determinations were then obtained directly from the measurement of the kinematical scale length in the mid-plane, 
or from the asymmetric drift relation. Here the analysis of the local kinematics of stars shows that this assumption is 
wrong in the Galaxy and introduces bias favouring much larger estimates of R p . 
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Recent analyses, using simpler distribution functions (Evans & Collett 1993 and Cuddcford & Binney 1992), one 
assuming R p = R a i the other a flat rotation curve, were not able to reproduce accurately the observed DF. Our results 
are concordant with the work by Fux & Martinet (1994) based on Jeans equations expanded by Cuddeford & Amendt 
(1991) with assumptions to obtain the closure of the moment equations. They deduce a short density scale length and 
they conclude, as we do here, that this result does not support the assumption <t z (R)/<tr(R) — const. 

The structural parameter values are obtained with a good accuracy. The small value of R p obtained with our 
analysis of the kinematics of neighbouring stars is consistent with most recent determinations of the density: Robin 
et al. (1992ab) and Ojha et al. (1996) based on star counts or Durand et al. (1996) based on COBE maps. We have 
shown that the hypothesis relating density and kinematic scale length R p = R a /2 is incorrect and that scale length 
determinations based on this assumption must be ruled out. 
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